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Abstract 

If electrons had an electric dipole moment (EDM) they would induce EDMs 
of atoms. The ratio of the atomic EDM to the electron EDM for a particular 
atom is called the enhancement factor, R. We calculate the enhancement 
factor for the francium and gold atoms, with the results 910 db ~ 5% for Fr 
and 260 ± ~ 15% for Au. The large values of these enhancement factors make 
these atoms attractive for electron EDM measurements, and hence the search 
for time-reversal invariance violation. 

PACS number(s): 11.30.Er, 32.10.Dk, 14.60.Cd 
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The existence of a non-zero electric dipole moment (EDM) of an atom, electron, or any 
quantum-mechanical system would imply that time-reversal invariance (T) and parity (P) 
are violated. To date, no non-zero EDM has been observed, though experimental limits on 
their magnitude have been obtained. The Standard Model predicts an upper limit on the 
electron EDM of the order of 10 _40 e cm ffl, while various alternative models predict values 
many orders of magnitude larger (for reviews of predicted values of the electron EDM see 
Refs. [0|§). Measurements of the electron EDM are worthwhile as in the future they may 
be able to distinguish between these models. Direct measurements of the electron EDM 
are difficult due to the electron's charge (see, e.g., [§|), so results for the electron EDM are 
obtained from atomic EDM measurements instead. The present limit on the electron EDM 
is \d e \ < 4 x 10 _27 e cm, from an experiment on the EDM of the Tl atom [|J. A summary of 
the results of atomic EDM measurements can be found in Table 6.2 of Ref. [[J. 

An atomic EDM can be induced by the presence of an electron EDM (as well as by other 
mechanisms, such as T- and P-violating electron-nucleon or nucleon-nucleon interactions, 
and nucleon EDMs). This allows us to obtain experimental results for the electron EDM 
by measuring the atomic EDM. The atomic EDM (d^) induced by the electron EDM (d e ) 
would be proportional to d e , and the ratio R = dj^jd e for a particular atom is known as the 
enhancement factor. As was first noted by Sandars ||, the atomic EDM of a heavy atom 
can be many times larger than the electron EDM; in fact, R is of the order of Z 3 a 2 times 
a relativistic factor (~ 3 in heavy atoms), where Z is the atomic number and a = 1/137 
(see, e.g., |7],|8|||). To convert experimental results for atomic EDMs to results for electron 
EDMs the value of the enhancement factor R is needed. Summaries of results for various 
enhancement factors can be found in Ref. [Q and Table 6.1 of Ref. ||. In this work we do 
accurate calculations of R for the Fr and Au atoms. 

The T- and P-odd interaction between the EDM of an electron and the electric field of 
the nucleus (— de'jo'S ■ E) results in an admixture of the ground state of the electron with 
excited states of opposite parity, according to perturbation theory. On one side of the atom 
the ground state and excited state wave functions will have the same sign, while on the 
other side they will have opposite signs. Therefore the total wave function will be larger on 
one side of the atom, hence the electron will be more likely to be there and so the atom will 
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have an EDM. An expression for this atomic EDM, and hence R, can be presented in the 
following form (see, e.g., Refs. 

fi= ^ = 2eERe (0N")("l(7.-l)S.E|0) (l) 

de n E — E n 

(for an atom with one valence electron), where |0) is the unperturbed ground state, {\n)} is 
the set of states with which it is mixed (including unbound, continuum states for which the 
sum should be replaced by an integral), — e is the charge on the electron, E is the electric field 
produced by the nucleus [= Zer/r 3 ; the main contribution to the second matrix element in 
(H|) comes from short distances, where the electric field of the nucleus is unscreened], z = r z , 
and 70 and S are the normal matrices of relativistic quantum mechanics. The ground state 
of the valence electron in Fr is 7s, while in Au it is 6s. The operator (70 — 1)S ■ E is a 
pseudoscalar and so it can only mix states having opposite parity and the same total angular 
momentum. Therefore the ground state is only mixed with pi/2 states. 
Equation (P can be rewritten as 
4Z« /o°°[/ s *(r)r/ np (r) + 9 :(r)rg np (r)}r 2 Q(r) dr J™ g* np (r)r~ 2 g s ( r yP(r) dr 

by using the following expression for the electron's relativistic wave function (see, e.g., ||): 
Vvzjm = (fnji(r)Qji m , -ig nj i(r)(cr ■ r/r)fi j7m ) T , where / and g are radial wave functions (f s 
refers to the ground state and f np the pi/ 2 excited states) and flji m is a spherical spinor 
(an eigenfunction of j 2 and j z ). Computer generated wave functions were used to calculate 
the integrals in Eq. (0) (where available, previously determined values were used). These 
wave functions were obtained using the relativistic Hartree-Fock method. The factors P(r) 
and Q(r) take into account the screening of the nuclear electric field by electrons and core 
polarization corrections (for the non-direct core polarization contribution these factors ac- 
tually become non-local operators). We also took into account correlation corrections to the 
wave functions. The many-body perturbation theory methods that we used are described in 
Refs. As a test, we also performed calculations using a semi-empirical method PJ 



that does not require computer calculations. The results were in good agreement with the 
numerical calculations. 

For Fr, we used the experimental value of the 7s-7pi/2 radial integral of r [the first integral 
in (fj)] that was determined in ||15|| : —5.238(10), in units of the Bohr radius (this compares 



well with the calculated value in fll4| : —5.241). For the 7s-8pi/2 radial integral of r we 
used the value calculated in |TJ| (we used the most complete many-body calculation value 
denoted by "Brueckner plus non-Brueckner" in Table IV of this work), with an estimated 
accuracy of 3%. The 7s-9pi/ 2 and 7s-10p\/ 2 radial integrals of r were calculated by us, as 
were the 7s-npi/ 2 values of the second integral in (Q) (all with an estimated accuracy of 
3%). We used the values of the 7s, 7pi/ 2 , 8pi/ 2 , and 9pi/ 2 energy levels listed in |16|], while 
we calculated the 10pi/2 energy level ourselves. We truncated the summation in (0) for the 
discrete states at the 10pi/ 2 state, as the remainder of this series is very small (it gives a 
contribution to the enhancement factor ~ 3, i.e., 0.3%). For the unbound, continuum states 
all integrals were calculated in the present work. For these we did not take into account 
screening, core polarization or correlation corrections, and so the errors for these integrals 
are larger (we estimate a 50% error for the whole continuum contribution), though this does 
not have an excessively large effect on the final error as the contribution of the continuum 
states to the enhancement factor is small (= 30). The final value of the enhancement factor 
for Fr is 910, with an estimated 5% error. 



For Au, we used the energy levels listed in [17]. Using the experimental result for the 



oscillator strength for the 6s-6pi/ 2 transition in [IS], we obtained the value of the 6s-6pi/ 2 
radial integral of r: —2.16(2). All of the other radial integrals were calculated in the present 
work, with an estimated accuracy of 10% for the discrete states. Once again, we estimate 
the error for the continuum contribution as 50% [the continuum contribution was again 
small (= 20)]. We truncated the summation over discrete states at the 8pi/ 2 state, with the 
remainder of the series giving a contribution ~ 1 (0.3%) to the enhancement factor. The 
final value of the enhancement factor for Au is 260, with an estimated 15% error. 

These results are in reasonable agreement with the previously determined estimates of 



the enhancement factors: ~ 1150 for Fr [|19|] and ~ 250 for Au 



4 



REFERENCES 

[1] M.E. Pospelov and LB. Khriplovich, Yad. Fiz. 53, 1030 (1991) [Sov. J. Nucl. Phys. 53, 
638 (1991)]. 



[2] 
[3] 
[4] 
[5] 

[6] 
[7] 



[9 
[10 

[11 
[12 

[13 

[14 

[15 

[16 

[17; 



W. Bernreuther and M. Suzuki, Rev. Mod. Phys. 63, 313 (1991). 
S.M. Barr, Int. J. Mod. Phys. A 8, 209 (1993). 

E.D. Commins, S.B. Ross, D. DeMille and B.C. Regan, Phys. Rev. A 50, 2960 (1994). 

LB. Khriplovich and S.K. Lamoreaux, CP Violation Without Strangeness (Springer- 
Verlag, Berlin, 1997). 

P.G.H. Sandars, Phys. Lett. 14, 194 (1965). 

V.V. Flambaum, Yad. Fiz. 24, 383 (1976) [Sov. J. Nucl. Phys. 24, 199 (1976)]. 

LB. Khriplovich, Parity Nonconservation in Atomic Phenomena (Gordon and Breach, 
Philadelphia, 1991). 

V.B. Berestetskh, E.M. Lifshitz, and L.P. Pitaevskh, Relativistic Quantum Theory 
(Pergamon Press, Oxford, 1971). 

V.A. Dzuba, V.V. Flambaum, PG. Silvestrov, and O.P Sushkov, J. Phys. B 20, 1399 
(1987). 

V.A. Dzuba, V.V. Flambaum, and O.P. Sushkov, Phys. Lett. A 140, 493 (1989). 

V.A. Dzuba, V.V. Flambaum, and O.P. Sushkov, Phys. Lett. A 141, 147 (1989). 

V.A. Dzuba, V.V. Flambaum, and O.P. Sushkov, Phys. Lett. A 142, 373 (1989). 

V.A. Dzuba, V.V. Flambaum, and O.P. Sushkov, Phys. Rev. A 51, 3454 (1995). 

J.E. Simsarian, L.A. Orozco, G.D. Sprouse, and W.Z. Zhao, Phys. Rev. A 57, 2448 
(1998). 



Reference [141 and references therein. 



C.E. Moore, Atomic Energy Levels, Natl. Bur. Stand. Ref. Data. Ser., Natl. Bur. Stand. 



(U.S.) Circ. No. 35 (U.S. GPO, Washington, D.C., 1971), Vol. III. 
[18] P. Hannaford, P.L. Larkins, and R.M. Lowe, J. Phys. B 14, 2321 (1981). 
[19] P.G.H. Sandars, Phys. Lett. 22, 290 (1966). 

[20] W.R. Johnson, D.S. Guo, M. Idrees, and J. Sapirstein, Phys. Rev. A 34, 1043 (1986). 



6 



